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1 Abstract 

2 Background 

3 Recently available whole genome sequences of three citrus species: one Citrus Clementina and two Citrus 

4 sinensis genomes have made it possible to understand the features of candidate disease resistance genes 

5 with nucleotide-binding sites (NBS) domain in Citrus and how NBS genes differ between hybrid and 

6 original Citrus species. 

7 Result 

8 We identified and re-annotated NBS genes from three citrus genomes and found similar numbers of NBS 

9 genes in those citrus genomes. Phylogenetic analysis of all citrus NBS genes across three genomes 

10 showed that there are three approximately evenly numbered groups: one group contains the Toll- 

1 1 Interleukin receptor (TIR) domain and two different groups that contain the Coiled Coil (CC) domain. 

12 Motif analysis confirmed that the two groups of CC-containing NBS genes are from different 

13 evolutionary origins. We partitioned NBS genes into clades using NBS domain sequence distances and 

14 found most clades include NBS genes from all three citrus genomes. This suggests that NBS genes in 

15 three citrus genomes may come from shared ancestral origins. We also mapped the re-sequenced reads of 

16 three pomelo and three Mandarin orange genomes onto the Citrus sinensis genome. We found that most 

17 NBS genes of the hybrid C sinensis genome have corresponding homologous genes in both pomelo and 

18 mandarin genome. The homologous NBS genes in pomelo and mandarin may explain why the NBS genes 

19 in their hybrid Citrus sinensis are similar to those in Citrus Clementina in this study. Furthermore, 

20 sequence variation amongst citrus NBS genes were shaped by multiple independent and shared 

21 accelerated mutation accumulation events among different groups of NBS genes and in different citrus 

22 genomes. 

23 Conclusion 
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1 Our comparative analyses yield valuable insight into the understanding of the structure, evolution and 

2 organization of NBS genes in Citrus genomes. There are significantly more NBS genes in Citrus 

3 genomes compared to other plant species. NBS genes in hybrid C. sinensis genomes are very similar to 

4 those in progenitor C. Clementina genome and they may be derived from possible common ancestral gene 

5 copies. Furthermore, our comprehensive analysis showed that there are three groups of plant NBS genes 

6 while CC-containing NBS genes can be divided into two groups. 

7 

8 Keywords: 

9 NBS genes, Comparative Genomics, Citrus Genomics 
10 

1 1 Background 

12 Disease resistance genes (R genes) are essential components of plant immune systems. Amongst five 

13 diverse classes of disease resistance genes [1], the largest class is those proteins that have nucleotide- 

14 binding site (NBS) domain and Leucine-Rich Repeat (LRR) domain. The NBS domain is evolutionarily 

15 conserved and is typically used to identify and characterize plant R genes. The LRR domains in R genes 

16 that mediate direct or indirect interactions [2, 3] with pathogen molecules are usually rapidly evolving to 

17 adapt to the change of pathogen ligands [4]. NBS domain-containing R genes were classified into two 

18 major types based on their domain structures in N terminus: proteins with a Toll-Interleukin receptor 

19 (TIR) domain and proteins with a Coiled Coil (CC) domain. While the TIR domain is usually well 

20 defined [5], the CC domain has higher variation and is less well characterized. 

21 Genome wide analyses of NBS genes in many genomes have showed that NBS R genes are 

22 diverse in number, structure and organization [5-12]. For example, grass genomes studied to date do not 

23 have TIR NBS genes; on the other hand, dicot genomes usually contain more TIR NBS genes than CC 

24 NBS genes [13]. Meanwhile, comparative analyses of NBS genes have facilitated the understanding of 

25 the organization, classification and evolution of NBS R genes. Comparison of R genes in multiple 
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1 genomes indicated that R genes are shaped by dynamic birth-and-death processes [14]. Comparative 

2 analyses of NBS genes in two Arabidopsis genomes (A. thaliana and A. lyrata) indicate that mating 

3 system shift from outcrossing to inbreeding has had a limited impact on the numbers of NBS genes, at 

4 least in short term[15]. Comparative analyses of NBS genes of diploid Phaseolus and tetraploid Glycine 

5 species concluded that whole genome duplication did not result in NBS R gene number increases. 

6 Recently, comparison of R genes in diverse grass genomes by Yang et al. showed that rapid evolution in 

7 R genes in Zea, Sorghum, and Brachypodium is associated with rice blast disease resistance [16]. 

8 Citrus species are amongst the most important fruit trees and have been cultivated for more than 

9 4000 years [17, 18]. Phylogenetic analyses using molecular markers showed that other cultivated citrus 

10 species, such as, sweet orange, grapefruit, and lemon, are derived from three original cultivated Citrus 

11 species: Citrus medica (citron), Citrus reticulata (Mandarin orange) and Citrus maxima (pomelo) [19, 20]. 

12 For example, the Citrus sinensis (sweet orange) is suggested to be the backcross hybrid of C. maxima 

13 (pomelo) and Citrus reticulata (Mandarin orange) [19-21]. Phylogenetic analyses of partial NBS genes 

14 of Poncirus trifoliata (trifoliate orange), Citrus reticulata (tangerine) and their Fl progeny showed that 

15 NBS genes of Poncirus trifoliata (trifoliate orange), Citrus reticulata formed genus-specific clades. 

16 Additionally, NBS genes of their Fl progeny had sister relationships to only one of the parents [22]. This 

17 suggests that NBS genes in crossing hybrid citrus species are also different from those in original citrus 

18 species. 

19 Recently available draft whole genome sequences of three citrus species: Citrus Clementina, 

20 Citrus sinensis from USA and Citrus sinensis from China [2 1 ], have made it possible to scan and identify 

21 all NBS genes in those genomes. In this study, we performed a genome-wide comparative analysis of 

22 NBS genes in three citrus genomes {Citrus Clementina, Citrus sinensis from China and Citrus sinensis 

23 from USA) to address the following questions: (1) What are the features of NBS genes in Citrus, such as 

24 numbers, physical locations, within-gene domain structures, and evolutionary dynamics of NBS genes?; 
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1 (2) Do NBS genes differ between hybrid and non-hybrid Citrus species?; and (3) Do Citrus NBS genes 

2 differ from NBS genes in other plant genomes? 

3 

4 Results 

5 Identification and classification of Citrus NBS Genes 

6 We searched the Citrus Clementina, Citrus sinensis USA and Citrus sinensis China genomes for genes 

7 containing the NBS domain using hmmsearch. Then, the NBS-containing genes were confirmed through 

8 homology searches against the Swiss protein database (see "Materials and Methods"). We identified 

9 similar numbers of NBS domain-containing genes amongst these three genomes. We found 618, 650 and 

10 508 NBS genes from Citrus Clementina, Citrus sinensis China and Citrus sinensis USA genomes 

11 respectively (Table 1). Among them, 413, 499 and 484 NBS genes were predicted in the original gene 

12 annotations and 205, 151 and 24 NBS genes in Citrus Clementina, Citrus sinensis China and Citrus 

13 sinensis USA respectively were newly predicted in this project (Table SI). 

14 NBS genes could be classified into different classes based on their domain structures [23]. We 

15 searched the Toll-Interleukin receptor (TIR) and Coiled Coil (CC) domains in the N-terminal region and 

16 the Leucine-Rich Repeat (LRR) domains in the C-terminal of NBS genes. We identified 117 NBS genes 

17 with CC-NBS-LRR (CNL) domains, 82 NBS genes with TIR-NBS-LRR (TNL) domains, 68 NBS genes 

18 with NBS-LRR (NL) and 351 others NBS genes without LRR domains (N, CN and TN) from Citrus 

19 Clementina (Table 1). We also identified 113 CNL, 77 TNL, 68 NL and 398 others NBS genes from 

20 Citrus sinensis China and 60 CNL, 30 TNL, 85 NL and 333 others NBS genes from Citrus sinensis USA 

21 (Table 1). In comparison to other genomes, there are many more Citrus NBS genes without the LRR 

22 domain. There are only 43.2%, 38.8% and 34.4% NBS genes with existing LRR domains in Citrus 

23 Clementina, Citrus sinensis China and Citrus sinensis USA, respectively, while 72% NBS genes in 

24 Arabidopsis [23] and 78% NBS genes mPopulus trichocarpa [24] have LRR domain. 
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1 The structures of TNL and CNL NBS genes are significantly different. The TNL NBS genes tend 

2 to have more introns than that of CNL NBS genes as previously found in Arabidopsis [23] and Populus 

3 [24]. The average numbers of introns in TNL NBS genes are 4.39, 4.70 and 4.47, while the numbers of 

4 introns in CNL NBS genes are 0.89, 1.15 and 1.38 in Citrus Clementina, Citrus sinensis China and Citrus 

5 sinensis USA respectively (Figure SI). The median numbers of introns of different types of NBS genes in 

6 the three Citrus genomes are similar: four in TNL NBS genes and one in CNL NBS genes. 

7 The NBS genes were unevenly distributed in the Citrus scaffolds/chromosomes. There were 217, 

8 158 and 110 NBS genes, totaling 78.3% of 618 NGS genes, distributed in scaffolds 5, 3, and 7 of C. 

9 Clementina, respectively. There were 125, 121 and 75 NBS genes distributed in chromosomes 1, 3 and 5 

10 of C. sinensis China. The majority (95 out of 107) of NBS genes with TIR domain in C. Clementina was 

11 distributed in scaffold 3. There were 33 and 51 (out of 102) NBS genes with TIR domain in C. sinensis 

12 China distributed on chromosome 5 and chromosome unknown, respectively. 

1 3 Phylogenetic and Clade Analysis of the Citrus NBS Genes 

14 We used the protein sequences of NBS domains, which are the most conserved part of NBS genes, to 

15 construct a phylogenetic tree of NBS genes. We only selected NBS domain sequences longer than 200 

16 amino acid residues and contain both P-loop and MHDV motifs. Finally, 442 Citrus Clementina, 393 

17 Citrus sinensis China, and 264 Citrus sinensis USA NBS domain sequences were used for phylogenetic 

18 analysis (Table SI). A Maximum Likelihood (ML) phylogenetic tree of these 1,099 NBS genes was then 

19 constructed using FastTree [25]. As shown in Figure 1, the un-rooted phylogenetic tree can be divided 

20 into three main groups. Most NBS genes containing the TIR domain are in one branch and most NBS 

21 genes containing a CC domain comprise the other two branches (Figure 1). Therefore, we denoted the 

22 three branches as TIR, CC1 and CC2. 

23 There were 452, 382 and 265 Citrus NBS genes in CC1, CC2 and TIR groups, respectively. 

24 Table S2 lists the group classification of 1,099 Citrus NBS genes. For each sub-group, we constructed a 

25 new ML phylogenetic tree rooted by Streptomyces coelicolor protein P25941 (Figure S2). The NBS 
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1 domains in CC groups, especially in the CC1 group, were relatively more diverged in sequence versus 

2 those in TIR group. The average Poisson corrected distances between sequence pairs within each group 

3 were 0.947, 0.827 and 0.663 for CC1, CC2 and TIR groups, respectively. We aligned the NBS domains 

4 sequences in each group using mafft with default parameters (—auto). The average percentages of 

5 identities for NBS domain sequences were 40%, 44% and 52% for CC1, CC2 and TIR groups, 

6 respectively. 

7 The number of LRR domains varied amongst the three classes of Citrus NBS genes. The number 

8 of LRR domains in the TIR group is significantly higher than the numbers of LRR domains in the CC 1 

9 and CC2 groups (Figure S3). The majority of CC1 NBS genes have only one LRR domain but the 

10 majority of TIR NBS genes have 2 LRR domains. The most frequent type of LRR domain in the CC1 

11 group is LRR 8 (PF13855), but LRR1 (PF00560) is the most frequent domain in both the CC2 group 

1 2 and the TIR group (Table S3). 

13 Previous studies showed that NBS genes from different species fall into separate phylogenetic 

14 groups [26]. In contrast, NBS genes in the three Citrus genomes we investigated are mixed amongst 

15 branches of phylogenetic tree as shown in the Figure 1. We partitioned the NBS gene tree into clades 

16 based on NBS domain sequence distance using PhyloPart with distance threshold of 0.025. This resulted 

17 in 1 14 clades (with more than two genes in each clade) and eight orphan genes (Table S2). The numbers 

18 of clades in each NBS gene group are similar. There are 39, 42, 32 clades in CC1, CC2 and TIR groups, 

19 respectively. We calculated average sequence similarities between NBS domains in each clade using 

20 alistat in SQUID (http://selab.janelia.org/software.html). The minimum average sequence similarities of 

21 NBS domains in each clade was 70%. The largest clade included 90 NBS genes. On average, there were 

22 9.56 genes per clade. There were seven clades with more than 40 genes. Among these clades, four clades 

23 were in CC1; two clades were in CC2 and one clade was in the TIR group. Most of the clades contain 

24 NBS genes from the three Citrus genomes (Figure 1, Table 2). There are 89 clades with more than three 

25 NBS genes and 87 of these clades contain members from each of the three Citrus genomes. Together, 
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1 these results imply that all of the NBS genes of the three Citrus genomes may have possible common 

2 ancestors. 

3 To understand the evolutionary dynamics maintaining functional constraint while maintaining so 

4 many gene family members, we tested for evidence of accelerated sequence evolution (positive selection) 

5 by searching for positively-selected sites within each clade. We detected positively selected sites in 53.5% 

6 (61) of clades. There were 18, 24 and 19 clades with sites under positive selection in CC1, CC2 and TIR 

7 groups, respectively. In total, we detected 541 positively selected sites. Consistent with previous reports 

8 [5], there are more positively selected sites in the C-terminal region of NBS genes (LRR domains) than 

9 those in the N-terminal and NBS domains (Table S4). Approximately 56.4% (305 out of 541) of the 

10 positively selected sites were located in the C-terminal region (LRR domain). While most clades have 

11 positively selected sites in the LRR domain, there were clades with a greater number of positively 

12 selected sites within the NBS domain. For example, there were 46 positively selected sites in Clade_1260 

13 that are located in NBS domain while there are only two positively selected sites in LRR domain (Figure 

14 S4). 

15 Mapping Re-sequencing Reads Showed the Conservation of NBS Genes in Citrus Genomes. 

16 To verify that Citrus sinensis (sweet orange) is a backcross hybrid of C. maxima (pomelo) and C. 

17 reticulata (Mandarin orange) [21], Xu et al. re-sequenced three C. maxima cultivars and three C 

18 reticulata cultivars and showed that SSR and SNP markers in the C. sinensis genome derive from the C. 

19 maxima and C. reticulata genomes at an approximately ratio of 1:3 [21]. Since our comparison of C. 

20 maxima and C. reticulata genomes showed that most of their NBS genes are very similar in sequence, we 

21 investigated the relationships amongst the NBS genes of C. sinensis (sweet orange), C. maxima (pomelo), 

22 and C. reticulata (Mandarin orange) to see if the evolutionary dynamics of NBS genes is informative. 

23 We mapped the re-sequenced reads of three Citrus maxima (pomelo) genomes and three C. 

24 reticulata (Mandarin orange) genomes to hybridized genome of Citrus sinensis (sweet orange) China 

25 using BWA with mem method [27] and calculated the coverage of each NBS gene in Citrus sinensis 
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1 China. We detected 524 of 650 Citrus sinensis China NBS genes with more than 50% coverage in all of 

2 the six resequenced genomes. Additionally, 249 of 650 Citrus sinensis China NBS genes with more than 

3 85% coverage in all six resequenced genomes and 592 (91%) of Citrus sinensis China NBS genes with 

4 coverage over 85% in at least one of the resequenced Citrus genomes. Using a cutoff of 15% coverage as 

5 missing calls over the entire length of the gene, only 25 (<4%) of Citrus sinensis China NBS genes seem 

6 to be deleted in at least one of the six resequenced genomes. No NBS genes were lost amongst all of the 

7 six resequenced genomes. 

8 Resequenced reads tended to map to exons rather than to introns (Table S5). For example, the 

9 entire coding region of NBS gene Cs9gl3310 was 3,693 bp in length and covered by all reads from the 

10 resequenced Citrus samples. However, the 16 kb introns were not consistently present in all of the 

11 resequenced reads. When we used exon presence instead of whole gene coverage, 403 out of 650 (62%) 

12 Citrus sinensis China NBS genes have more than 85% coverage in all six resequenced Citrus genomes 

13 and 645 (99%) Citrus sinensis China NBS genes have more than 85% coverage in at least one of the 

14 resequenced Citrus genomes. Most Citrus sinensis China NBS genes have high exon coverage on both 

15 the C maxima (pomelo) and C. reticulata (Mandarin orange) genomes (Figure S5). There were only three 

16 Citrus sinensis China NBS genes, Cslg02140.1, Cs6g02120.1 and Cs7g02220.1, with low (<40%) exon 

17 coverage in all three C. maxima genomes, but high coverage was achieved in at least two of the three C. 

18 reticulata genomes (Figure S5). Our mapping results showed that most NBS genes of C. sinensis (sweet 

1 9 orange) have a corresponding copy in the resequenced C. maxima (pomelo) and C reticulata (Mandarin 

20 orange) genomes. 

2 1 Motif Patterns in Citrus NBS Genes 

22 To further examine the gene structures of three Citrus NBS gene groups, we searched the motifs in the 

23 sequences of each group using MEME. Figure 2 lists the top 20 motifs identified by MEME from the 

24 CC 1 , CC2 and TIR Citrus NBS gene groups. 
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1 The motifs in the N-terminal of CC 1 and CC2 groups showed little similarity to each other. The 

2 MEME found two motifs from CC1 group. One of them, EEQQQMRRLNQVQGWLSRVEA, was 

3 present in 356 of 452 CC1 NBS genes. The other motif, GSQEIDKLCLGGYC SKNCKS S YKFGKKVA, 

4 was present in 189 of 452 CC1 NBS genes. Both motifs have high level of sequence similarities with the 

5 CC NBS genes of Arabidopsis (Figure S6 (A)). The MEME found three motifs from CC2 group. Two of 

6 the motifs, KKLTNMLEMIKAVLDDAEEKQ and AVKLWLGKLKDAAYDVEDVLDEFQTEALR 

7 were identified in 325 and 340 of 382 CC2 NBS genes, respectively. The motif 

8 AVKLWLGKLKDAAYDVEDVLDEFQTEALR has high levels of similarity with the motif identified 

9 from the CC NBS genes of Oryza sativa japonica (Japonica rice) [28] (Figure S6 (B)). The different 

10 motifs in CC1 and CC2 NBS genes implied that they may be from different evolutionary origins. 

1 1 For the TIR group, the MEME identified five motifs in the N-terminal. They can be found in 76% 

12 to 90% of 265 TIR Citrus NBS genes. The first four motifs YDVFLSFRGEDTRDNFTSHLY, 

13 AIEASAISVIIFSEGYASSRWCLDELVKI, GQIVIPVFYRVDPSDVRKQTG, and 

14 ENPEKVQKWRDALKEA were very similar to the TIR1-4 motifs identified from TIR NBS genes in 

15 Arabidopsis [23] and Populus trichocarpa [24] (Figure S6 (C)). 

16 The MEME algorithm identified eight motifs from NBS domains of CC2 and TIR groups, which 

17 are similar to the motif structures of NBS domains in Arabidopsis [23] and Populus trichocarpa NBS 

18 genes [24]. MEME results also showed that the motif structure of NBS domain of CC1 groups is slightly 

19 different from those of CC2 and TIR groups. Ten motifs were identified from NBS domains of CC1 

20 group with two extra RNBS motifs (RNBS-F, RNBS-G) between GLPL and MHDV motifs. Five of 

21 eight motifs: P-loop, Kinase-2, RNBS-B, GLPL and MHDV, showed high levels of sequence similarity 

22 amongst the three NBS groups, which were also similar to those motifs from Arabidopsis [23] and 

23 Populus trichocarpa NBS genes [24] . The MHDV motif in Citrus was often slightly modified to 

24 MHDL, as found in Arabidopsis and Populus previously. Meanwhile, the motifs RNBS-A, RNBS-C, and 

25 RNBS-D were quite dissimilar to each other amongst three groups. 
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1 The LRR domains in the C-terminal of NBS genes usually have high sequence diversity as they 

2 play a role in recognizing pathogen virulence proteins. The MEME algorithm identified five, six, and four 

3 motifs from LRR domains of CC1, CC2 and TIR NBS genes. Most of the motifs contained LxxL repeats. 

4 Motifs from different groups were highly variable in sequence, which implied that the three groups of 

5 NBS genes have different levels of sequence, and thus functional, constraint and play different roles in the 

6 Citrus immune system. Some motifs had repeated several times in the same gene. For examples, motif 

7 APNLKSLEVSSCxxMEEIISV was found in 367 of 452 CC1 group with an average 4.8 motifs per gene; 

8 motif IKTLPESVCELYNLQTLDLEGCRRLKKLP was found in 332 of 382 CC2 group with an average 

9 2.8 times per gene and motifs CKRLKSLPSSLCKLKSLGxLNLSGCSNLE and 

10 GNISELFLDGTAIEELPS SIE were found in 208 and 209 of 265 TIR group with 2.5 and 2.1 times per 

1 1 gene, respectively. 

12 Analysis of NBS Gene Clusters in Citrus Genomes 

13 The majority of Citrus NBS genes were physically clustered in genome (Table 3). 525 of 618 NBS genes 

14 in C. Clementina were found in 108 clusters and 500 of 650 NBS genes in C. sinensis China were found in 

15 126 clusters (Table S2). Although the assembly of C. sinensis USA is more fragmentary, there still are 

16 207 of 508 NBS genes present in 72 clusters. The largest number of gene clusters in C. Clementina, C 

17 sinensis China and C. sinensis USA contain 55, 37 and 13 NBS genes, respectively. 

18 Most clusters contain NBS genes from the same group. In C. Clementina, there were 38 clusters 

19 with NBS genes of CC1 group, 40 clusters with NBS genes of CC3 group and 23 clusters with NBS 

20 genes of TIR group. Only seven out of 108 clusters contain the NBS genes from two or three groups. In 

21 C. sinensis China, there were 41 clusters with NBS genes of CC1 group, 44 clusters with NBS genes of 

22 CC2 group and 26 clusters with NBS genes of TIR group. There were 15 clusters containing NBS genes 

23 from two or three groups in C. sinensis China. In C. sinensis USA, there were 27 clusters with NBS genes 

24 of CC1 group, 25 clusters with NBS genes of CC2 group and 15 clusters with NBS genes of TIR group. 
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1 Only five clusters contain NBS genes from different groups in C. sinensis USA. The lower number of 

2 clusters in C. sinensis USA may due to its fragmental assembly. 

3 The sequences of NBS genes within clusters are much more similar to each other than those 

4 between clusters (T-test p value < 2.2e-16). The mean identities between genes within and between 

5 clusters were 0.628 and 0.393 respectively (Figure 3). Furthermore, NBS genes in the same cluster tend to 

6 be in the same strand, which indicates that the NBS genes in the clusters are due to tandem duplication. 

7 We detected 204 pairs of tandem duplications in C. Clementina and 217 pairs of tandem duplications in C. 

8 sinensis China using MCScanX [29]. The numbers of tandem duplications within CC1 and CC2 groups 

9 are much greater than that within the TIR group. Among 204 pairs in C. Clementina, 90 and 85 pairs were 

10 present in the CC1 and CC2 groups and only 29 pairs from TIR group. Among 214 tandem gene pairs in 

11 C. sinensis China, 88 and 87 pairs were present in the CC1 and CC2 groups and only 42 pairs in the TIR 

12 group. The fewer tandem duplications of NBS genes in the TIR group may be the reason that there are 

1 3 fewer clusters of TIR groups in Citrus genomes. 

14 We identified 254 and 246 gene conversion events from 76 NBS gene clusters in C. Clementina 

15 and 75 NBS gene clusters in C. sinensis China, respectively (Table 4). The gene conversion events in C. 

1 6 sinensis USA is much less due to fragmental assembly. It is interesting that most of the conversion events 

17 (483 out of 520) were identified from the relatively small clusters with less than 10 NBS genes. Among 

18 these conversions, 119 events located in N-terminal, 178 in NBS domains and 223 in C-terminal, which 

1 9 indicating that there was no significantly bias in the location of conversion. Most gene conversions were 

20 between genes from the same group. We identified 101, 184 and 229 conversion events from NBS genes 

21 in TIR, CC1 and CC2 groups respectively and only four conversion events were identified between NBS 

22 genes of TIR and CC2 groups. While we could identify similar total numbers of conversion events in C. 

23 Clementina and C. sinensis China, we found almost double such events in C. sinensis China TIR clusters 

24 versus in C. Clementina TIR clusters. Meanwhile, we also found many more conversion events in C. 

25 Clementina CC2 clusters than that in C. sinensis China CC2 clusters (Table 4). 
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1 Analysis of NBS Orthologs 

2 We identified 719 Citrus NBS gene pairs of orthologs amongst Citrus Clementina, C. sinensis China and 

3 C. sinensis USA. 270 orthologs were shared between C Clementina and C. sinensis China (Figure S7); 

4 227 orthologous gene pairs were shared between C. Clementina and C. sinensis USA and 222 orthologous 

5 gene pairs were shared between C. sinensis China and C. sinensis USA. The percentages of identities 

6 between orthologous genes range from 53.8% to 100% (mean of 93.6%, median of 96.8%). The 

7 percentages of identities between orthologous genes from TIR group (TNL and TN, 88.42 ± 10.04) were 

8 significantly (T-test: P<2.6e-6) lower than those from CC groups (CNL and CN, 92.92 ± 7.03). The 

9 nonsynonymous divergence (dN) values of the 719 orthologs ranged from 1.46e-6 to 0.31 (mean of 0.03, 

10 median of 0.08) and the synonymous divergence (dS) values ranged from 1.14e-5 to 0.75 (mean of 0.05, 

1 1 median of 0.019) (Figure 4). The dN/dS ratios ranged from 0.001 to 50 with a median of 0.62. Similar to 

12 Arabidopsis, the dN and dS of orthologs from TIR group (TNL and TN) were higher than those from CC 

13 groups (CNL and CN). The dN values of orthologs from TIR group ranged from 0.001 to 0.24 (mean of 

14 0.038, median of 0.017) while those of orthologs from CC group ranged from 1.4e-6 to 0.23 (mean of 

15 0.025, median of 0.007). The dS of orthologs from TIR group ranged from 3.27e-5 to 0.76 (mean of 

16 0.078, median of 0.029) while those of orthologs from CC group ranged from le-5 to 49 (mean of 0.04, 

17 median of 0.016). However, the dN/dS ratios of orthologs from TIR group were relatively lower than that 

18 of orthologs from CC group (Figure 4(A)). 

19 The dN and dS rates as well as dN/dS rate ratios of ortholog residues in clusters are generally 

20 greater than those of orthologous singletons (Figure 4(B)). The median dN/dS ratio of 148 orthologous 

21 singletons was 0.55 while that of 361 orthologs in clusters was 0.63. There were 124 orthologs with 

22 dN/dS ratios above 1. These orthologs could undergo positive selection pressure. Most of the positive 

23 selected orthologs (87.9%, 109 out of 124) belonged to CC groups (49 orthologs from CC1 group and 50 

24 from CC2 group). 
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1 We also detected 38 NBS gene syntenic blocks between C. Clementina and C. sinensis China 

2 using MCScanX [29] (Figure S8). On average, there are 1 1 genes per block. The biggest block contains 

3 45 NBS genes and was located in scaffold_5 of C. Clementina and chromosome 3 of C. sinensis China. 

4 There were 80, 77 and 71 syntenic NBS genes in scaffold_5 of C. clementinal 'chromosome 5 of C. 

5 sinensis China, scaffold_3 of C. clementina/chromosome 5 of C. sinensis China and scaffold_7 of C. 

6 clementina/chromosome 1 of C. sinensis China. Furthermore, there were 68, 66 and 8 syntenic NBS 

7 genes in the unknown chromosome of C. sinensis China corresponding to NBS genes in scaffold three, 

8 five and seven of C. Clementina, respectively. 

9 Analysis of Conserved NBS Gene Clusters 

10 We identified 118 pairs of conserved NBS gene clusters between C. Clementina and C. sinensis China (Figure 

11 S9). There were 19 NBS gene clusters completely conserved in C. Clementina and C. sinensis China. For 

12 example, all 14 NBS genes in cluster CL225 of C. sinensis China have orthologs in the cluster CL168 C. 

13 Clementina which contains nine NBS genes and vice versa. Some clusters (in both C. Clementina and C. 

14 sinensis China) have several corresponding conserved clusters. This may due to either the genome arrangement 

15 or the incomplete assembly of genomes. Furthermore, there were 19 clusters of C. Clementina with no 

16 conserved clusters in C. sinensis China and 18 clusters of C. sinensis China with no conserved clusters in 

17 C. Clementina. 

18 The conserved clusters provided additional insights into NBS gene evolution within and between 

19 Citrus genomes. For example, cluster CL142 (9 NBS genes) in C. Clementina and cluster CL282 (7 NBS 

20 genes) in C. sinensis China were highly conserved. The phylogenetic tree of these 16 genes suggested 

21 division into two subgroups (depicted in blue and orange color in Figure 5 (A)). This division indicates 

22 two ancestral genes for this conserved cluster. Using the phylogenetic tree as a framework, we 

23 reconstructed the evolutionary history of these two clusters. There were several tandem duplication events 

24 in the evolutionary history of the conserved clusters, and an extra tandem duplication was observed in C. 

25 Clementina after it separated from C. sinensis (Figure 5 (B)). Two NBS genes were lost in C. Clementina 
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1 and one NBS gene was lost in C. sinensis. Furthermore, there was a recombination event within C. 

2 Clementina. 

3 Mutations and Transposons in Citrus NBS Genes 

4 Plant NBS genes are continuously evolving. Sequence variation and structural constraints are shaped by 

5 gene birth-and-death processes [30, 31]. Besides possible interallelic recombination and gene conversion, 

6 gene mutations and transposable elements appear to play important roles in NBS gene evolution. 

7 We compared the DNA sequences of NBS genes from different Citrus genomes to identify 

8 mutations. Approximately half of Citrus NBS genes have mutations maintained amongst corresponding 

9 orthologs. When a mutation took place in an exon and resulted in stop-codon gaining or frame-shift, the 

10 target gene often became a pseudogene. For example, Cslgl8610.1_cc_32, Cslgl8610.1 and 

11 orange l.lg003 3 67m are orthologs of C. Clementina, C. sinensis China and C. sinensis USA. There were 2 

12 stop-codon gaining mutations in Cslgl8610.1, which means that Cslgl8610.1 might become a pseudo 

13 (nonfunctional) gene after gaining these mutations. 

14 LTR retrotransposons are widespread in eukaryotic genomes, especially plant genomes. We 

15 predicted 19014, 6296 and 1479 LTR retrotransponsons with typical LTR characters in the draft genomes 

16 of C. Clementina, C. sinensis China and C. sinensis USA, respectively. Then, we filtered out LTR 

1 7 retrotransponsons with low similarity with known TE proteins using BLASTX with e-value greater than 

18 le-5. Finally, 4920, 3726 and 1240 LTR retrotransponsons remained in C. Clementina, C. sinensis China 

19 and C. sinensis USA, respectively. We identified 33, 32 and 4 NBS genes that were inserted with LTR 

20 retrotransposons in C. Clementina, C. sinensis China and C. sinensis USA respectively (Table S6). Most 

21 of these genes will likely become pseudogenes due to these insertions. For an instance, 

22 orangel.lg043039m, orangel.lg043039m_cc_l 16 and orange l.lg043 03 9m_csc_ 123 were orthologs of 

23 C. sinensis USA, C. Clementina and C. sinensis China respectively (Figure S10 (A)). The results of gene 

24 structure analysis showed that structure in orange l.lg043 03 9m in C. sinensis USA seems relatively well 

25 maintained. But there is an about 10 kb fragment of LTR-retrotransposons in the corresponding 
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1 homologous gene in C. Clementina and C. sinensis China (Figure S10 (B)). We identified both 

2 orange l.lg043 03 9m_cc_l 16 and orangel.lg043039m_csc_123 as pseudogenes using PseudoPipe [32]. 

3 Experimental Validation of One NBS Gene in Citrus Genomes 

4 We validated the orthologs of a conserved NB gene, Cslg09350.1, in a wide range of Citrus species. 

5 Cslg09350.1 conserved in the 3 sequenced Citrus genomes. It is a CNL NBS gene and has 5 exons with 

6 about 4 kb in length. We designed the primers targeting about 3.5 kb fragment and could amplify a 3.5 kb 

7 fragment using PCR in different Citrus species, including C. sinensis (sweet orange), C. Clementina 

8 (clementine), C. japonica (Kumquat), C. sinensis Navelina (Navel orange), C maxima (pomelo), C 

9 aurantiifolia (lime), and C. aurantium (sour orange). The orthologs from C. sinensis and C. Clementina 

10 were identical as expected. There are a few mutations in the orthologs of C. japonica, C. sinensis 

1 1 Navelina (Navel orange), and C maxima, but the orthologs of C. aurantiifolia and C aurantium should 

12 be undergoing a pseudogenization process. There was a deletion plus several mutations in the NBS 

13 domain in C. aurantiifolia. There was an eight-base deletion in the LRR domain in C. aurantiifolia and C. 

14 aurantium. Most Citrus aurantiifolia mutations were shared with C. aurantium (Figure Sll), and these 

15 mutations were not shared with other species. The results suggest that in C. aurantiifolia and C. 

16 aurantium, this gene is likely derived from the same common ancestor and was inherited as a pseudogene. 

17 

18 Discussion 

19 Possible Common Ancestor of NBS Genes in Hybrid Citrus sinensis and Original Citrus Clementina 

20 Genome 

21 After carefully reannotating the Citrus genome sequences, we found similar numbers of NBS genes in 

22 Citrus Clementina and Citrus sinensis China. There are slightly fewer NBS genes in Citrus sinensis USA, 

23 possibly due to the more fragmental assembly of this genome. In phylogenetic tree using NBS domains, 

24 the NBS genes from three different Citrus genomes are mixed together. After partitioning NBS genes on 

25 the tree into clades, we found that 97.7% of the clades containing more than three genes had members 
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1 from each of three Citrus genomes. This pattern suggests that these three Citrus genomes have similar 

2 NBS genes derived from common ancestors. Because Citrus sinensis is the hybrid of C. reticulata 

3 (Mandarin orange) and C maxima (pomelo), it should be heterozygous and some of NBS genes in Citrus 

4 sinensis are expected have different genetic distances from NBS genes in Citrus Clementine. This would 

5 support previous observations of NBS genes in Fl progeny of Poncirus trifoliata (trifoliate orange) and 

6 Citrus reticulata (tangerine). However, this is clearly not the case. 

7 Furthermore, we mapped the resequenced reads of three C maxima (pomelo) genomes and three 

8 C. reticulata (Mandarin orange) genomes onto the genome of Citrus sinensis China. In this case, 62% of 

9 C sinensis China NBS genes have a copy present in all six resequenced genomes and 99% of C. sinensis 

10 China NBS genes have a copy in at least one of the resequenced genomes. The mapping results confirmed 

11 that a significant percentage of NBS genes of hybrid C. sinensis genomes have corresponding 

12 homologous genes in both the C. maxima and the C. reticulata genomes. Because the reference genome 

13 sequence of C. maxima is not yet available, the total number of NBS genes in C. maxima genomes is still 

14 not known. However, we can at least conclude from the mapping of resequenced genomes that the C. 

15 maxima genome has homologous copies of NBS genes in C. reticulata and C. sinensis genomes. The 

16 homologous NBS genes in C. maxima and C. reticulata may be the reason that NBS genes in their hybrid 

17 Citrus sinensis are similar to those in C. reticulata in this study. 

1 8 Three Groups of Citrus NBS Genes 

19 We identified 442, 393 and 264 genes with full length NBS domains from Citrus Clementina, Citrus 

20 sinensis China and Citrus sinensis USA reference genomes, respectively. There are also many genes with 

21 short NBS domains in three Citrus genomes. The Citrus NBS genes can be divided into three groups 

22 according to the phylogenetic tree of NBS domains: two of them contain CC domain and the other group 

23 contains TIR domain. The number of CC NBS genes is three times of the number of TIR NBS genes. In 

24 most of the TIR NBS genes, we can find the LRR domains defined in Pfam database [33]. We only can 

25 identify LRR domains in small part of CC NBS genes using Pfam LRR domain definition. However, we 
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1 can find the LxxLs repeats in most of the CC NBS genes as shown in motifs from MEME. This implied 

2 that there may be other types of LRR domain in Citrus CC NBS genes. Our motif analyses also showed 

3 that motifs of TIR domains in Citrus CC2 and TIR NBS genes are similar to those of TIR domains in 

4 Arabidopsis [23] and Populus trichocarpa [24] TIR NBS genes. Furthermore, the motifs of Citrus CC 

5 domains in CC1 NBS genes are similar to those motifs of CC domains in Arabidopsis CC NBS genes and 

6 the motifs of Citrus CC domains in Citrus CC2 NBS genes are similar to those motifs of CC domains in 

7 japonica rice CC NBS genes [28]. The different structure of motifs in NBS domain and CC domain 

8 between Citrus CC1 NBS genes and Citrus CC2 NBS genes implied that they are from different 

9 evolutionary origin. 

10 To further confirm the three groups of NBS genes, we identified the NBS genes from 

1 1 Arabidopsis, Populus, Oryza sativa and grape. We used the same criteria to select the NBS domains from 

12 NBS genes of these genomes. Finally, we selected 152, 216, 209, 126 NBS domains from Arabidopsis, 

13 Populus, Oryza sativa and grape, respectively. Then, together with 442 NBS domains selected from 

14 Citrus Clementina, we constructed a phylogenetic tree of those 1145 NBS domain sequences from five 

15 genomes. As shown in Figure 6, the un-rooted phylogenetic tree was divided into three main branches. 

16 The Populus, grape and Citrus genomes have significant amount of NBS genes in all three branches. 

17 However, the NBS genes of Oryza sativa dominated in CC2 branch and most of Arabidopsis NBS genes 

18 located in CC1 and TIR branches. Our study showed that the NBS genes can be divided into three major 

19 groups as the NBS genes with the CC domains are separated into two groups. The three groups of NBS 

20 genes underwent divergent evolution in different genomes. Further comparison of NBS genes of more 

21 genomes may help to understand the evolution of NBS genes and will help elucidate how plants maintain 

22 and adapt their defense system against pathogens. 

23 Highly Clustering of Citrus NBS Genes 

24 The Citrus NBS genes are highly clustered in the genome. 84.9% of NBS genes in Citrus Clementina and 

25 76.9% of NBS genes in Citrus sinensis China were found in clusters. Previous studies showed 76% of 
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1 rice [34], 64% of A. lyrata and 71% A thaliana [15], 83.2% of grapevine and 67.5% of poplar [11] NBS 

2 genes are found in clusters. The percentage of Citrus NBS genes in clusters is in the high level comparing 

3 to other genomes. The average number of genes per clusters is 4.86 in Citrus Clementina and 3.97 in 

4 Citrus sinensis China. These numbers are similar to those in other genomes [11]. Furthermore, most 

5 Citrus NBS genes in the same cluster belong to the same phylogenetic group, suggesting that tandem 

6 duplication is the primary mechanism for the expansion of NBS genes in the Citrus genus. 

7 Molecular Evolution of Citrus NBS Genes 

8 Similar to NBS genes in other genomes, Citrus NBS genes are highly dynamic and are shaped by several 

9 evolutionary processes leading to several differences amongst NBS genes, including domain presence and 

10 mutation constraints and genome organization. Our results revealed multiple molecular evolution events 

1 1 amongst Citrus NBS genes including gene duplications, gene conversions, mutation constraint changes, 

12 recombination and transposable element insertions. Likely these events support a birth-and-death process 

13 leading to the origins of new NBS genes as well as malfunction and loss of other NBS genes [14]. We 

14 found more than 200 tandem duplications in both C. Clementina and C. sinensis China genomes alone. 

15 We also found that NBS genes become pseudogenes following original frame-shift mutations leading to 

1 6 mutation accumulations plus disruption of gene constraint leading to loss of function through transposable 

1 7 element insertions. 

18 Most molecular evolutionary events occurred within NBS gene groups, suggesting that gene 

19 birth-death processes have been occurring since divergence from a common ancestral gene copy. 

20 Interestingly, the molecular evolution processes occurred differently amongst NBS gene groups and 

2 1 differently within each Citrus genome. For example, the number of gene conversion events in C sinensis 

22 China TIR clusters is almost double that in C. Clementina TIR clusters, while there are many more 

23 conversion events in C. Clementina CC2 clusters than in C. sinensis China CC2 clusters (Table 4). 

24 Numbers of tandem duplications in CC1 and CC2 groups are much greater than that from the TIR group. 

25 
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1 Conclusions 

2 Our comparative analyses yield valuable insight into the understanding of the structure, evolution and 

3 organization of NBS genes in Citrus genomes. There are significantly more NBS genes in Citrus 

4 genomes compared to other plant species. Citrus NBS genes are structurally highly clustered. NBS genes 

5 in hybrid C. sinensis genomes are very similar to those in progenitor C. Clementina genomes and their 

6 NBS genes may be derived from possible common ancestral gene copies. Furthermore, our 

7 comprehensive analysis also showed that there are three groups of plant NBS genes while NBS genes 

8 containing CC domains can be divided into two groups. 

9 

10 Materials and Methods 

1 1 Sequences Used 

12 We downloaded the draft genome sequences and the original gene annotations of Citrus Clementina 

13 (clementine) and Citrus sinensis USA (sweet orange from USA) from the Citrus Genome Database 

14 (http://www.citrusgenomedb.org/) and those of Citrus sinensis China (Chinese sweet orange [21]) from 

15 the Citrus sinensis annotation project (http://citrus.hzau.edu.cn/orange/). The sizes of assembled genomes 

16 of Citrus Clementina, Citrus sinensis China and Citrus sinensis USA are 301, 328 and 319 million 

17 basepairs, respectively. The C. Clementina genomes assembled into nine major scaffolds and 95.8% of the 

18 sequences were assigned to those nine scaffolds. About 72.9% of the C. sinensis China genome assigned 

19 to nine chromosomes. The genome of C. sinensis USA only assembled to 12,574 scaffolds and the N50 of 

20 scaffolds is 250 kb. The original gene annotations of Citrus Clementina, Citrus sinensis China and Citrus 

21 sinensis USA have 24,533, 29,385 and 25,397 genes respectively. We also downloaded the resequence 

22 data of three Citrus Clementina (clementine) and three Citrus maxima (pomelo) genomes from the Citrus 

23 sinensis annotation project (http://citrus.hzau.edu.cn/orange/). 

24 Identification of NBS Genes in Citrus Genomes 
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1 We first screened the original predicted citrus open reading frames (ORFs) using hmmsearch [35] with 

2 the hidden Markov models (HMM) of Pfam [36] for NBS domain presence (NB-ARC, PF00931) using 

3 an e-value cut-off of 0.1 for the hmmsearch. Then, the proteins selected by the HMM were searched 

4 against the Swiss protein database [37] to confirm the annotation using BLASTP [38]. Only the proteins 

5 that have a significant match (e-value < 1E-5 in BLASTP search) with the NBS proteins or resistant 

6 proteins in the Swiss protein database were identified as NBS-containing proteins. To recover possible 

7 NBS genes that may be missed in the original gene annotations, we mapped the identified NBS genes to 

8 the draft genome using TBLASTN. The matched sequences with e-value < 1E-5 were then predicted 

9 using Genewise [39]. The new genes predicted by Genewise were also confirmed using NBS domain 

10 HMM screening and by BLASTP searching through the Swiss protein database. 

1 1 Identification of Orthologous NBS genes in C. Clementina and C. sinensis 

12 Orthologous NBS genes of C. Clementina and C. sinensis were identified using the reciprocal best blast 

13 method [40]. We used NBS genes of C. Clementina as query sequences to search against the NBS genes 

14 of C. sinensis and vice versa. Protein pairs with reciprocal best hits of e-value < 1E-20 were defined as 

15 orthologs. 

16 To calculate the rates of nonsynonymous, synonymous and their rate ratio (dN, dS and dN/dS) of 

1 7 orthologous pairs, we first aligned orthologous protein sequence pairs using rnafft, and then converted the 

18 protein alignments to codon-based alignments using PAL2NAL [41]. We calculated the dN, dS and 

19 dN/dS ratios using the codeml program in PAML version 4.7 [42]. 

20 Phylogenetic Analysis 

21 We constructed the phylogenetic tree of NBS genes of the three citrus genomes using only the conserved 

22 NBS domain. Only 1,099 NBS domain sequences that have both the P-loop and the MHDV motifs and 

23 were longer than 70% (200 amino acids) of the full-length NBS domain were included. Next, NBS 

24 domain sequences were aligned in mafft [43] using an auto alignment model and a best fit maximum 

25 likelihood phylogenetic tree of NBS genes was constructed using FastTree [25] with default parameters 
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1 (JTT+CAT). The resultant best fit phylogenetic tree was divided into 3 main groups based upon clade 

2 support. For each group, we constructed a new ML tree with Streptomyces coelicolor protein P25941 as 

3 the outgroup using FastTree [25]. The average identity of each group was calculated using alistat 

4 implemented in SQUID (http://selab.janelia.org/software.html). We further partitioned the NBS gene tree 

5 into clades using the depth-first phylogeny partition method in PhyloPart [44] with distance threshold 

6 0.025. This clustered the NBS genes into 1 14 clades, which contain more than 2 genes, and 8 orphan 

7 genes. 

8 Domain and Motif Annotation 

9 The Toll-Interleukin receptor (TIR) domains in Citrus NBS-containing proteins were identified using 

10 hmmsearch [35] with the HMM model of Pfam domain PF01582 and a 0.1 e-value cut-off. The Leucine- 

1 1 Rich Repeat (LRR) domains in Citrus NBS containing proteins were identified using the HMM models of 

12 Pfam LRR domains with e-value cut-off of 0.1. As long as there is a significant hit to one of LRR 

13 domain models (e-value < 0.1), it was define as an LRR-containing protein. We used MARCOIL [45] 

14 with a threshold probability of 90 and COILS [46] with a threshold of 0.9 to search for Coiled Coil (CC) 

15 domains in the N-terminal region of Citrus NBS-containing protein. We considered a protein as a CC- 

16 containing protein if either MARCOIL or COILS reported a CC domain in it. 

17 We identified 20 motifs amongst the NBS genes in each of the three main phylogenetic groups 

18 separately using MEME SUITE [47]. The motif width was set to between 6 and 50 for MEME. Then, we 

1 9 searched the motif structure of all genes in each group using MAST with default parameters (-ev 1 0 -mt 

20 0.0001). 

21 Pseudogene Identification 

22 We identified possible pseudogenes using PseudoPipe [32] with default parameters (-e 0.1). The 

23 PseudoPipe algorithm identifies pseudo genes by integrating sequence similarity, intron-exon structure, 

24 plus presence of stop codons and frame-shifts. We used all Citrus NBS genes to search Citrus genomes 

25 for potential NGS pseudo genes. 
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1 Transposon Identification 

2 All long terminal repeat (LTR) retrotransponsons in each citrus genome were identified using LTR finder 

3 [48] with default parameters (-0 3 -t 1 -e 1 -m 2 -u -2). Then, we used a script program to match the 

4 location of the LTR transposons to the NBS-LRR genes in each citrus genome. 

5 NBS Gene Synteny Identification 

6 Gene synteny was identified and defined using MCScanX [29] with default parameters (-A -u 5000). First, 

7 NBS genes between two genomes were aligned using BLASTP and the matches with E-value < le-5 were 

8 sorted according to their chromosome positions. Synteny scores were then calculated for each block based 

9 upon gene position. Two genes were considered in the same block if there were fewer than 25 genes 

10 separating them. MCScanX reported blocks with at least 5 collinear gene pairs. 

1 1 Gene Cluster Analysis 

12 We grouped the NBS genes in each Citrus genome into the same cluster if the genome location between 

13 two genes was within 200 kb. We also identified the conserved gene clusters between the C. Clementina 

14 and C. sinensis. If all genes in a cluster of C. Clementina have orthologs in the corresponding cluster of C. 

15 sinensis and vice versa, then these two clusters were called completely conserved clusters. If only part of 

1 6 genes in the clusters have orthologs, we called these two clusters partially conserved. 

17 Gene Conversion Detection 

18 We first aligned the sequences of NBS genes in the same cluster using rnafft [43]. Then, we used 

19 GENECONV [49] version 1.81a with default settings (N=l 0,000) to detect gene conversions. 

20 GENECONV identifies gene conversions by finding identical fragments between pairs of sequences in a 

21 nucleotide alignment. A global P value < 0.05 was used to assess the statistical significance of the 

22 observed conversions. GENECONV requires at least three sequences for analyses in order to account for 

23 shared ancestral states. Thus, we only detected conversions in clusters containing three or more genes. 

24 Tests for Sites under Positive Selection 
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1 The amino acid sequences of NBS genes from the same clade were aligned with mafft [43]. Then we 

2 converted the protein alignments to codon-based alignment using PAL2NAL [41]. The positively selected 

3 sites were statistically identified using the Bayesian approach implemented in codeml within PAML [42]. 

4 We also further examined sites in the co > 1 class with >90% posterior probability. 

5 Mapping of Re-sequencing Data of Citrus Clementina and Citrus maxima 

6 We mapped the raw reads of each re-sequenced sample to the draft genome of Citrus sinensis China using 

7 BWA [27] with default parameters (-k 19 -d 100 -A 1 -B 4 -O 6). Then, the mapped reads of NBS genes 

8 regions were extracted using BEDtools [50]. Two types of coverage of each NBS gene in Citrus sinensis 

9 were calculated. One coverage type divides the length of mapped sequences by the whole gene sequence 

10 length and the other type divides the length of mapped sequences by the length of exon sequence only. 

1 1 Citrus DNA Extraction and PCR Amplication 

12 Citrus leaf samples were collected from the six Citrus plants in USHRL's (USDA Horticultural Research 

13 Laboratory, Fort Pierce, Florida): Citrus sinensis (sweet, Navel orange), Citrus aurantium (Karum jamir, 

14 sour orange), Citrus reticulata (Mandarin orange), Citrus Clementina (Clementine), Citrus aurantiifolia 

15 (sweet lime), Citrus japonica (Yuzu, kumquat), and Citrus maxima (pomelo). Total DNA was extracted 

16 from leaf midribs following the DNeasy® Plant Mini Kit standard protocol (Qiagen Inc., Valencia, CA), 

17 followed by DNA quantity and quality evaluation with Nanodrop. We chose the NBS gene, Cslg09350, 

18 which was conserved in C. Clementina and C. sinensis for validating the conservation of NBS gene 

19 among different Citrus genomes. Primers used in this study were designed using Oligo 7.23 (Molecular 

20 Biology Insights, Inc., Cascade, CO, USA). High Fidelity Platinum® Taq DNA Polymerase (Invitrogen, 

21 Carlsbad, CA, USA) was used to amplify the NBS-LRR genes from Citrus DNA. For PCR, 20 uL 

22 reactions using standard conditions provided by the manufacturer for High Fidelity Platinum® Taq DNA 

23 Polymerase. PCR was performed using an initial denaturation at 95°C for 3 minutes, 35 cycles of 94°C 

24 for 20 seconds, 50-52°C for 20 seconds (specified by different primer sets) and 68°C for 3 minutes, follow 

25 by final extension at 68°C for 10 minutes in a CI 000™ Thermal Cycler (Bio-Rad, Hercules, CA). The 
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1 cloning and sequencing analysis of amplified PCR products were conducted as previously 

2 described [51] 
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8 

9 

10 

11 Figure Legends: 

12 Figure 1. Maximum likelihood phylogenetic tree of Citrus NBS-LRR genes constructed from 

13 multiple sequences alignment of NBS domain. 

14 There are three groups in the phylogenetic tree: CC1, CC2 and TIR. Clades were classified using 

15 PhyloPart as shown in alternating color. The outer circle shows species with Citrus Clementina (Cc) in 

1 6 red, Citrus sinensis China (CsCN) in green and Citrus sinensis USA (CsUSA) in blue. 

17 

1 8 Figure 2. Architecture of NBS gene motifs in Citrus. 

19 Top 20 motifs identified by MEME for each group of citrus NBS genes are listed. The motifs underlined 

20 in green were repeated in C-terminal (LRR domain). 

21 

22 Figure 3. Percentage identity distribution of NBS-encoding genes in C. Clementina (Cc) and C. 

23 sinensis China (Csc). 

24 Green lines indicate pairwise identity distribution of inter-cluster NBS-encoding genes. Orange lines 

25 indicate intro-cluster NBS-encoding genes. 

26 
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1 Figure 4. Divergence of orthologous NBS-encoding genes amongst Citrus species. 

2 A) Divergence of different classes. B. Divergence of singletons and cluster members. 
3 

4 Figure 5. Duplication histories of one of the conserved NBS-encoding genes cluster in C. Clementina 

5 and C. sinensis. A. 

6 Phylogenetic tree of the NBS-encoding genes; B. Orthologs in the conserved cluster of C. Clementina and 

7 C. sinensis. Black connectors indicate tandem duplications and red dotted connectors indicate gene loss. 

8 A gene recombination event occurred in the ancestor of C. Clementina, as indicated with an orange 

9 connector. 

10 

1 1 Figure 6. Phylogenetic analysis of NBS genes of Citrus Clementina, Arabidopsis thaliana, Populus 

12 trichocarpa, Oryza sativa and Vitis (grape). 

13 Species are indicated by the color of the outer circle and branches. 
14 

15 Tables: 

16 Table 1. Classification of citrus NBS genes 



Class 


Citrus Clementina 


Citrus sinensis China 


Citrus sinensis USA 


CNL 


117 


113 


60 


TNL 


82 


77 


30 


NL 


68 


62 


85 


Others 


351 


398 


333 


Total 


618 


650 


508 



17 

18 The CNL class of NBS genes contains three domains: Coiled Coil (CC), nucleotide-binding sites (NBS) 

19 and Leucine-Rich Repeat (LRR). The TNL class of NBS genes contains three domains: Toll-Interleukin 

20 receptor (TIR), nucleotide-binding sites (NBS) and Leucine-Rich Repeat (LRR). The NL class of NBS 

21 genes contains two domains: nucleotide-binding sites (NBS) and Leucine-Rich Repeat (LRR). The other 

22 NBS genes contain no LRR domain. 

23 
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Table 2. Type of citrus NBS gene clades 

NBS genes from citrus genomes Number of Clade 

Citrus Clementina, Citrus sinensis China, Citrus sinensis USA 87 

Citrus clementina, Citrus sinensis China 13 

Citrus clementina, Citrus sinensis USA 7 

Citrus sinensis China, Citrus sinensis USA 5 

Citrus clementina, 1 

Citrus sinensis China, 0 

Citrus sinensis USA 1 



Table 3. Citrus NBS genes in clusters 



Citrus species 



Total NBS Number of NBS genes 



genes 



in cluster 



Number of Percent of NBS genes in 
clusters cluster (%) 



Citrus clementina 618 
Citrus sinensis China 650 
Citrus sinensis USA 508 



525 
500 
207 



108 
126 
72 



84.9 
76.9 
40.7 



Table 4. Gene conversion events found in C. clementina and C. sinensis China 







Number of clusters 


Number of 




Mean size of 


Citrus species 


Cluster type 


(with conversion 
events) 


conversion 
events 


Affected genes 


conversion 
tracts 




TIR 


27 


35 


32 


111 


Citrus 












clementina 


CC1 


43 


86 


67 


137 




CC2 


46 


133 


83 


146 




TIR 


36 


69 


42 


137 


Citrus sinensis 












China 


CC1 


51 


89 


48 


162 




CC2 


57 


92 


75 


116 
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1 Support Information 

2 Supplemental Figures 

3 Figure SI. Average intron number of CNL and TNL. Cc: C. Clementina, CsCN: C. sinensis China and 

4 CsUSA: C. sinensis USA. 

5 

6 Figure S2. Phylogenetic trees of three groups of citrus NBS gene with P25941 as outgroup. A, CC1 

7 group; B, CC2 group and C, TIR group. 

8 

9 Figure S3. Structure of domains of citrus NBS genes. The red rectangles indicated the NBS domains, the 

10 green ellipses indicated CC domains, the left pointing pentagrams indicated TIR domains and the left 

1 1 pointing triangle indicated LRR domains. 

12 

13 Figure S4. Positive selection sites of Clade_1260. The grey box indicated the NBS domain. 
14 

15 Figure S5. The exon coverage of NBS-encoding genes in the resequences of citrus genome (3 mandarin 

16 and 3 pomelo). 

17 

18 Figure S6. Multiple alignments of motifs from CC1 domain (A), CC2 domain (B) and TIR domain (C). 
19 

20 Figure S7. Orthologs of citrus NBS-encoding genes of Cc and CsCN and their genome locations. The 

21 outer blue circles indicated 9 largest scaffolds of Cc and orange ones indicated 9 chromosomes plus an 

22 un-located pseudo-chromosome. The NBS-encoding genes were arranged due to their locations on the 

23 chromosomes. The NBS genes in the TIR group were indicated with orange links, CC1 in blue and CC2 
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1 in light blue ones. The NBS-encoding orthologs of C. Clementina and C. sinensis were indicated with the 

2 link lines. 

3 

4 Figure S8. Dot plotting of syntenic orthologs of the NBS genes between C. Clementina (Cc) and C. 

5 sinensis China (CsCN) produced by MCScanX. The numbers in the grids indicated the number of 

6 syntenic ortholog pairs in the corresponding scaffolds. The scaffolds of C. Clementina (Cc) were arranged 

7 in x-axis and that of C. sinensis China (CsCN) were arranged in y-axis. 

8 

9 Figure S9. Diagram of conserved cluster between C. Clementina and C. sinensis. The nodes with orange 

1 0 color indicated the clusters of C. Clementina and the blue ones indicated clusters of C. sinensis and size of 

1 1 the nodes were relative to cluster size. The numbers on the edges indicated the number of orthologs could 

12 be found from the other species. 

13 

14 Figure S10. Pseudogenization of NBS-encoding genes of citrus due to mutation (A) and retrotransposon 

15 insertion (B). A: Sequences alignment of Cslgl8610.1, Cslgl8610.1_cc_32 and orange l.lg003 3 67m. 

16 The bases indicated the variation bases in each citrus species. The star in red color represented stop-codon 

17 gaining mutation and may result in pseudogene of Cslgl8610.1. B: LTR retrotransposon insertion in 

18 NBS-encoding genes orange l.lg043 03 9m_cc_l 16 and orange 1 . 1 g043 03 9m_csc_ 123 from Cc and CsCN 

19 respectively, and the multiple sequences alignment of the orthologs. LTR, Long terminal repeats; TSR, 

20 Target site repeats. 

21 

22 Figure Sll. Validation and phylogenetic analysis of NBS-LRR gene, Cslg09350.1, in different citrus 

23 species. A) multiple sequences alignment, the blue bars indicate the exons and the thin lines indicate the 

24 introns. The vertical lines in each alignment blocks indicate mutations comparing with the reference 
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1 sequences from C. sinensis. B) neighbor-joining phylogenetic tree of the orthologs of Cslg09350.1 in 

2 different citrus species and the bootstrapping values displayed on the branches. 

3 

4 Supplemental Tables 

5 Table SI. The NBS genes in three Citrus genomes 

6 Table S2. Classification, cluster and clade of Citrus NBS genes 

7 Table S3. LRR domain distribution of citrus NBS gene 

8 Table S4. Positive selection sites distribution 

9 Table S5. The exome coverage of NBS genes from the 3 Clementine and 3 pummelo resequencing 

1 0 genomes 

1 1 Table S6. Citrus NBS gene interapted with LTR transposons 
12 
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cc PAC 20814924 

cscCs1g13460.1 

cc PAC 20800606 

csu PAC 18094653 

cc PAC 20799641 

cscCs1g13270.1 

csu PAC 18097646 

esc Cs1g13380.1 

cc PAC 20801 600 

csu PAC 18097659 

cscorange1.1g038398m esc 103 

csu PAC 18134431 

cc PAC 20800321 

csu PAC 18134409 

cscorange1.1t01365.1 

cc PAC 20800907 

csu PAC 18105550 

cc PAC 20800869 

csu PAC 18116737 

esc clementineO.9 esc 56 

cc PAC 20800054 

csu PAC 18103212 

cc PAC 20800452 

csu PAC 18115065 

cc PAC 20801841 

cc PAC 20800033 

- esc Cs1g11430.1 

csu PAC 18115094 

cscorange1.1t01362.1 

cc PAC 20799565 

csu PAC 18103214 

cscCs1g11320.1 

csu PAC 18116030 

csu PAC 18115098 

cc PAC 20801943 

cscorange1.1t01435.1 

csu PAC 18115091 

csu PAC 18115092 

cscCs1g11390.1 

cc PAC 20801660 

cc PAC 20800484 

cc PAC 20800357 

cscCs1g11410.1 

csu PAC 18115119 

cc PAC 20801 002 

cc PAC 20800156 

csu PAC 18096063 

cc PAC 20812573 

esc Cs5g28645.2 

cscCs7g14510.1 

csu PAC 18129731 

csu PAC 18113378 

csu PAC 18113377 

cc PAC 20804331 

cc PAC 20803137 

cc PAC 20804269 

esc Cs7g14540.1 

csu PAC 18129698 

cc PAC 20801970 

csu PAC 18116732 

cc PAC 20799746 

cc PAC 20800575 

csu PAC 18115953 

cc PAC 20802277 

cc PAC 20802136 

cc PAC 20799757 

cscCs1g12190.1 

cc PAC 20800232 

csu PAC 18120966 

esc clementine0.9 esc 19 

cc PAC 20801961 

cc PAC 20799402 

cc PAC 20799564 

cc PAC 20802212 

cscCs1g12140.1 

cscorange1.1t01359.1 

csu PAC 18133381 

cc PAC 20800173 

esc orange1.1g036236m esc 95 
—csu PAC 18134301 
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cc PAC 20799363 

cc PAC 20800299 

cc PAC 20799768 

cc PAC 20800692 

cc PAC 20802305 

cc PAC 20802101 

csu PAC 18093006 

cc PAC 20800638 

cscCs1g12050.1 

cscorange1.1t01356.1 

cc PAC 20800410 

csu PAC 18134419 

cscorange1.1t04526.1 

cscCs1g12210.1 
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cc PAC 20801425 

esc clementine0.9 esc 17 

cc PAC 20794482 
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csu Cs1g07180.1 csu 8 
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cscCs1g07220.1 
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esc Cs8g07700.2 
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esc Cs8g07870.1 
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cc PAC 20817794 
cc PAC 20815783 
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csu PAC 18092884 - - 

cc PAC 20785660 

csu PAC 18115687 -- 
csu PAC 18126616 - - 
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cc Ciclev10003121m cc 2 

cc orange1.1g047163m cc 134 - 

esc Cs3g08140.1 

cscorange1.1t03126.1 

cc PAC 20786071 

cc PAC 20787804 

esc Cs3g08040.1 

cc orange1.1t03712.1 cc 152 -■ 

cc PAC 20789293 

esc Cs3g07940.2 

cc PAC 20786188 

cc orange1.1t03137.1 cc 146 -■ 

cscorange1.1t03707.1 

cc Ciclevl 0000338m cc 1 

csu PAC 18124435 

csu PAC 18124479 

cc PAC 20810732 

csu PAC 18124413 

esc Cs5g22030.4 

csu PAC 18113661 

cc PAC 20785477 

csu PAC 18124146 

csu PAC 18113653 

cscCs3g13160.2 

cc PAC 20785475 

cc PAC 20788943 

cscorange1.1t04538.1 

cc PAC 20786815 

esc Cs3g12720.1 

cc PAC 20787045 

esc orange 1 . 1 104537. 1 

csu PAC 18098033 

cc PAC 20786767 

cscCs3g12710.1 

csu PAC 18117062 

esc Cs3g12850.1 

cc PAC 20787446 

cc PAC 20788237 

cc PAC 20787088 

cc orange1.1t04537.1 cc 188 -■ 

esc Cs3g12760.1 

cc PAC 20818774 

ccCs3g12810.1 cc 56 

esc Cs3g 12810.1 

cc PAC 20785454 

cc PAC 20788609 

cscorange1.1t04078.1 

cc PAC 20789147 

esc clementine0.9 esc 43 

cc PAC 20786759 

esc orange 1 . 1 104094. 1 

cscorange1.1t04089.1 

cscorange1.1t04090.1 

cc PAC 20789280 

cc PAC 20786949 

esc clementine0.9 esc 58 

csu PAC 18098041 

esc Cs3g1 1190.1 

cc PAC 20788476 

- csu Ciclev10000331m csu 1 

cscCs3g11010.1 

cc PAC 20787724 

cc PAC 20786776 

cc PAC 20788790 

cc PAC 20788777 

cc PAC 20787177 

csu PAC 18115479 

cc PAC 20788413 

cc PAC 20787062 

cscCs3g11080.1 

cc PAC 20786766 

cscorange1.1g043056m esc 124 
esc Cs3g11030.1 
cscorange1.1t04079.1 
cc PAC 20788284 

- csu PAC 18108787 

- - cc PAC 20797441 

cscCs3g11210.1 
cc PAC 20786937 
cc PAC 20789277 
cscorange1.1t04084.1 
cc PAC 20787871 

- cc PAC 20788493 
csu PAC 18128264 
csu PAC 18128260 
esc Cs5g26170.1 
cc Cs5g26170.1 CC64 

-cc PAC 20812056 
-cscCs5g26190.1 
esc Cs8g09140.1 
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csu PAC 1813481 
esc Cs3g13250.1 
cc orange1.1t04930.1 cc 205 

- cc PAC 20788084 

esc Cs3g06650.1 
esc clementine0.9 esc 
cc PAC 20786990 
csu PAC 18103403 
cc PAC 20799321 
cc PAC 20786012 
cc PAC 20799325 
cc PAC 20787885 

- - cc PAC 20785870 
esc Cs3g06590.1 
cc Cs3g06650.1 cc 51 
cc PAC 20789463 

- cc PAC 20786523 

cscorange1.1t04831.1 

esc Cs3g13430.1 

cc PAC 20788009 

csu PAC 18111830 

esc Cs6g02120.1 

cc Cs6g02130.1 cc 67 

cc PAC 20796583 

cc PAC 20798378 

csu PAC 18113412 

esc Cs6g02100.2 

esc Cs6g02130.1 

- cc PAC 20797980 

cc PAC 20786190 

cscCs3g13200.1 

cscCs3g13210.3 

cc PAC 20787622 

cc PAC 20788708 

csu PAC 18103407 

csu PAC 18129839 

esc Cs3g06640.1 

cc Cs3g06490.1 cc 50 

cc orange1.1g003996m cc 93 

cscorange1.1t04833.1 

cscorange1.1t04933.1 

esc Cs3g06490.1 

cscorange1.1t04653.1 

cc PAC 20789465 

esc Cs3g06500.1 

cc PAC 20785487 

esc Cs3g06600.1 

cc PAC 2078581 1 

cscorange1.1t04830.1 -- 
cc PAC 20785472 
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ccCs3g06480.1 cc49 

cc PAC 207871 13 

cc PAC 20786065 

cc PAC 20788855 

cc orangel .1t04830.1 cc 202 
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Figure S2-B 
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P25941 

cc Cs3g05870.1 cc 44 - 

— csu PAC 1 81 35661 - - 

cc PAC 20795145 

csu PAC 1810991 1 

csu PAC 18114029 

esc Cs4g08 1 1 0.2 

cc PAC 20795174 

cc PAC 20796050 

cc PAC 20794844 

cc Cic[ev10010550m cc 7 

csu PAC 18133666 

esc orangel .1t03742.2 

cscCs1g17140.1 

cc PAC 20799749 

cc PAC 20813877 

csu PAC 18124524 

esc Cs8g09010.1 

cc PAC 20817766 

cc PAC 20816114 

csu PAC 18103803 

esc orangel. 1g002750m esc 84 

csu PAC 18124147 

cc PAC 20818998 

cc PAC 20815684 

csu PAC 18124119 

cscCs2g30580.1 

cc PAC 20800344 

cscCs1g24520.1 

csu PAC 18098638 

cc PAC 20807901 

cc PAC 20808995 

esc Cs8g12970.1 

esc Cs8g12930.1 

cc PAC 20808304 

cc PAC 20806949 

cc PAC 20806024 

cc PAC 20810036 

esc Cs1g18610.1 

ccCs1g18610.1 cc 32 

csu PAC 18103672 

cc Cs1g18380.1 cc 31 

esc Cs1g18370.1 

cc Cs1g18380.1 cc 30 

cscCs1g18380.2 

csu PAC 18133031 

cc Cs1g16990.1 cc 27 

esc Cs1g16990.1 

csu PAC 18101360 

csu PAC 18111753 

esc orangel. 1g037574m esc 101 

cc PAC 20799505 

csu PAC 18132959 

cc Cs1g17000.1 cc28 

cscCs1g17000.1 

cscCs1g17130.1 

csu PAC 18133159 

ccCs1g17130.1 cc 29 

csu PAC 18120763 

cc PAC 20800831 

ccCs1g15550.1 cc 26 

csu PAC 18120804 

cscCs1g15550.1 

csu PAC 18131396 

esc Cs1g09330.1 

ccCs1g09330.1 cc 24 

csu PAC 18117024 

esc Cs1g01890.1 

cc Cs1g01890.1 cc 15 

cc Cs1g01880.1 cc 14 

esc Cs1g01880.1 

csu Cs1g09350.1 csu 12 

cscCs1g09350.1 

cc Cs1g09350.1 cc 25 

csu PAC 18103079 

cscCs1g01180.1 

cc Cs1g01180.1 cc 13 

cc orangel. 1g042336m cc 113 - 

csu PAC 18099873 

esc clementine0.9 esc 31 

ccCs9g13120.1 cc 72 

csu PAC 18093623 

esc Cs9g 131 20. 1 

esc orangel. 1g038751m esc 107 

csu PAC 18134368 

cc orangel. 1g038751m cc 106 - 

csu PAC 18093759 

cc orangel. 1g038751m cc 105 - 

esc Cs1g09470.1 

csu PAC 18103075 

cc Cs1g01160.1 cc 12 

cscCs1g01160.1 

csu PAC 18113135 

csu PAC 18103101 

cc Cs1g01140.1 cc 11 

cscCs1g01140.1 

esc clementine0.9 esc 40 

csu PAC 18134367 

cc orangel. 1g002044m cc 79 - 
cc orangel. 1g043039m cc 116 - 

csu PAC 18134371 

esc orangel. 1g043039m esc 123 

esc Cs3g 1 3320. 1 

cc Cs3g13320.1 cc 57 

esc orangel. H02518.1 

esc Cs3g13390.1 

cc Cs3g13740.1 cc 60 

esc Cs3g13740.1 

cc Cs3g13740.1 cc 61 

csu PAC 18124354 

ccCs3g13710.1 cc 58 

esc Cs3g 13710.1 

csu PAC 18104953 

esc Cs3g13690.1 

ccCs3g13730.1 cc 59 

cscCs3g13730.1 

csu PAC 18094768 

esc Cs3g13340.1 

csu PAC 18127287 

csu PAC 18115706 

cc orangel. 1104513.1 cc 186 

csu PAC 18119372 

esc orangel. H04656.1 

esc orangel. H04515.1 

cc orangel. 1g002265m cc 82 - 

cc PAC 20786683 

csu PAC 18129247 

csu PAC 18119374 

esc clementine0.9 esc 29 

cc Cs3g04020.1 cc 39 

cc orangel. 1104656.1 cc 199 

csu PAC 18098231 

cc orangel. 1104656.1 cc 198 

esc clementine0.9 esc 28 

csu PAC 18105557 

esc Cs3g04020.1 

esc Cs3g03980.1 

cc orangel. 1104659.1 cc200 -■ 

esc orangel .1t04659.1 

esc clementine0.9 esc 33 

esc clementine0.9 esc 34 

cc orangel. 1g002071m cc 80 - 

csu PAC 18110658 

esc orangel. H01942.1 

esc orangel. 1101918.1 

cc orangel. 1101918.1 cc 142 - ■ 
cc orangel. 1101926.1 cc 143 - ■ 

esc orangel. H01926.1 

csu PAC 18116825 

csu PAC 18128569 

esc orange1.1t05161.1 

esc orangel .1t03739.1 

cc orangel. 1103734.1 cc 160 - ■ 

csu PAC 18128568 

cc orangel. 1103732.1 cc 155 - ■ 
cc orangel. 1103732.1 cc 154 
cc orangel. 1103732.1 cc 153 

csu PAC 18098063 

esc orangel .1103732. 1 

cc orangel. 1103732.1 cc 157 

cc orangel. 1101898.1 cc 141 - ■ 

esc orangel .1105262.1 

esc orangel. H03734.1 

esc clementine0.9 esc 7 

cc orangel. 1103732.1 cc 156 
cc orangel. 1103732.1 cc 159 - ■ 
cc orangel. 1103732.1 cc158 -■ 
esc orangel .1101 898.1 

csu PAC 18098236 

cscCs5g21975.1 

cscCs5g28770.1 

csu PAC 1 8096069 

cc PAC 20812381 

cc PAC 20788596 

esc Cs3g16680.1 

csu PAC 18107536 

csu PAC 18115546 

csu PAC 18100952 

esc Cs1g08080.2 

cc PAC 20795881 

cc PAC 20790459 

csu PAC 18116027 

esc Cs9g 1 3200. 1 

csu PAC 18123812 

esc Cs9g13250.1 

csu PAC 18128683 

cscCs9g13350.1 

cc PAC 20790057 

csu PAC 18108555 

esc orangel. 1g042863m esc 122 

cc PAC 20792539 

csu PAC 18123291 

cc PAC 20789492 

cc PAC 20789725 

esc clementine0.9 esc 12 

esc Cs9g13290.1 

esc Cs9g13300.1 

esc orangel. 1g037039m esc 98 

cc PAC 20796217 

cc PAC 20790303 

csu PAC 18103274 

cscCs9g13310.1 

cc PAC 20802397 

cscCs7g16030.2 

csu PAC 18110905 

esc Cs7g02220.1 

csu PAC 18120915 

cc PAC 20802693 

csu PAC 18100379 

cc PAC 20788086 

esc Cs6g14690.1 

cc PAC 20797957 

csu PAC 18118131 

csu PAC 18130652 

esc orangel. 1g043647m esc 127 

cc PAC 20794674 

csu PAC 18130926 

esc orangel. 1g048135m esc 139 

cc PAC 20796205 

csu PAC 18125888 

esc Cs4g07730.2 

esc orangel. 1t03487.1 

csu PAC 18105257 

cc PAC 20789997 

esc orangel .H03486.1 

cc PAC 20790721 

csu PAC 18117736 

esc orangel. 1g039990m esc 113 

csu PAC 18108759 

esc orangel .1103481 .1 

csu PAC 18117737 

esc clementine0.9 esc 2 

cc PAC 20807819 

esc orangel .1102478.1 

csu Ciclev10023525m csu 4 

cc PAC 20811240 

cc PAC 20808815 

esc orangel .1102481 .1 

esc orangel. H02477.1 

csu PAC 18092976 

esc orangel .1t02498.1 

cc PAC 20810956 

csu orangel .1102498.1 csu 29 - 

cscCs5g20560.1 

cc PAC 20807498 

cc PAC 2081 1577 

csu PAC 18115717 

cc PAC 20806592 

esc orangel. 1g001281m esc 82 

csu PAC 18122236 

esc orangel .1102495.1 

cc PAC 20807247 

cc PAC 20808198 

esc Cs5g20390.1 

esc Cs5g20430.1 

esc Cs5g20470.1 

cc PAC 20808565 

csu PAC 18104567 

esc Cs5g20500.1 

csu PAC 18104557 

cscCs5g20580.1 

cc PAC 20810018 

esc orangel .1t02493.1 

csu PAC 18100702 

cc PAC 20808822 

esc clementineO.9 esc 6 

cc PAC 20808683 

csu PAC 18119824 

esc orangel. 1104582.1 

esc Cs5g21 120.1 

csu PAC 18121255 

cscCs5g21040.1 

cc PAC 20810353 

csu Cs5g21040.1 csu 16 

cc PAC 20809006 

cc PAC 20808287 

esc clementine0.9 esc 53 

cc PAC 20808452 

cc PAC 20809129 

esc orangel .H04584.1 

esc clementine0.9 esc 16 

csu PAC 18103200 

csu PAC 18110894 

csu PAC 18113113 

cc orangel. 1g045012m cc 122 - 

esc clementine0.9 esc 50 

esc clementine0.9 esc 37 

csu PAC 18129815 

csu PAC 18134878 

esc Cs6g01070.1 

cc orangel. 1g037340m cc 101 - 

csu PAC 18122381 

esc Cs6g01080.1 

csu PAC 18103199 

esc Cs6g01040.1 

csu PAC 18123958 

csu PAC 18120192 

esc Cs1g05730.1 

cc Cs1g05730.1 cc 17 

cc PAC 20810431 

esc orangel .1t03534.1 

csu PAC 18124438 

cc Cs6g01020.1 cc 66 

esc Cs6g05250.1 

esc Cs6g05240.1 

cscCs6g01020.1 

csu PAC 18124427 

esc orangel. 1g038448m esc 105 

csu PAC 18113799 

cc Cs5g29510.1 cc 65 

;Cs5g29510.1 

;Cs5g22730.1 

j PAC 18122176 

PAC 20808955 

j PAC 18122178 

cscCs5g22710.1 

cc PAC 20805880 

cc PAC 20785531 

esc orangel. H01159.1 

csu PAC 18116086 

csu PAC 18116075 

esc clementine0.9 esc 1 

cscCs3g01180.1 

csu PAC 18123113 

cc orangel. 1103874.1 cc 165 - ■ 

esc orange1.1t03647.1 

csu PAC 18137817 

cc orangel. 1103647.1 cc 150 - ■ 
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